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The index distribution in the cross section of a multimode fiber has an 
important influence on the modal group velocities and, hence, on the fiber 
impulse response. In this paper we derive a method for the evaluation of 
arbitrary circular symmetric index profdes. In particular, we compute the 
impulse response of a fiber with a ring-shaped parabolic index profile 
which exhibits useful equalizing properties. The pulse spread is found to 
be nearly one order of magnitude smaller than that of a fiber with an equal, 
but abrupt, index decline from core to cladding. 

I. INTRODUCTION 

Multimode operation of optical fibers relaxes the fabrication toler- 
ances, allows the use of incoherent sources, and can alleviate handling 
and splicing problems. Modal (group) delay differences are nearly 
equalized 12 if the core index decreases as the square of the fiber radius 
from a maximum at the axis (Fig. 1). A distribution of this kind is 
realized in the Selfoc* fiber, which was indeed reported to have very 
low values of differential mode delay. 3,4 

Since then, the question has been raised whether there are other 
index profiles which have similar equalizing effects, but are otherwise 
perhaps more amenable to certain fabrication techniques or have 
advantages with respect to splicing or bending. Although the latter 
part of this question is difficult to answer at this time, it is certainly 
possible to identify at least one profile that has quite effective equaliz- 
ing properties. Imagine a slab with a square-law index distribution in 
transverse direction. The group velocities of all its modes are known to 
be nearly equal. 1 It is then plausible to expect that these properties 
are approximately preserved if the slab is warped in a way which 
results in a tube with the cross-sectional index distribution shown in 
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Fig. 1 — Concentric parabolic index profile. 

Fig. 2. Assume that a cladding material of lower refractive index fills 
the bore and surrounds the tube to the outside, so that a fiber is formed 
which guides modes within a tube-like structure with parabolic index 
distribution. 

The purpose of this paper is to identify the modes of this structure, 
calculate their group velocities, and predict the impulse response to 
be expected when all modes propagate uncoupled and with equal 
power. To do this, we employ the WKB description 5 in a form which 
ignores the anomalies of dielectric waveguide modes near cutoff, as- 
suming that few of all the propagating modes are close to this condition. 
For the sake of simplicity, we also restrict the following computations 
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Fig. 2 — Cross section Uumigh a ring-shaped parabolic index profile. Maximum 
index no along a circle of radius Re. Cladding index rio(l — A). 
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to small index variations, so that all propagation directions can be 
assumed paraxial to the waveguide axis and the corresponding ap- 
proximations apply. 

II. A CHARACTERISTIC EQUATION FOR CIRCULAR SYMMETRIC INDEX 
DISTRIBUTIONS 

Let us adopt a cylindrical coordinate .system (r, 0. z) and assume 
that the refractive index n is a function of r only. We define a local 
wave number 

k(r) = 2vn{r)/\ ) (1) 

where X is the wavelength in free space. Because of the circular sym- 
metry, we can separate the general wave equation and solve for and 
z. In doing so, we define an axial propagation constant /3 and describe 
the azimuthal periodicity by an azimuthal mode number v. The re- 
maining partial differential equation for the radial field dependence 
E(r) has then the form 



Following the usual WKB approach, 5 we substitute 

E(r) = e«<'>, (3) 

ignore the second derivative dhi'dr-, and, by solving for du/dr, we 
obtain the solution 



Given /3 and v, we can find two radii, Ri and R*, at which the root in 

(4) vanishes (Fig. 3). These radii define a ring-shaped region within 
which eq. (4) has an imaginary part causing the field E to be a periodic 
function. Outside of the region, E decreases or increases ^periodically. 

As in the 2-dimensional case, 5 decreasing (or evanescent) field charac- 
teristics outside are obtained if the total phase inside the region is 

[ ' Vfc 2 (r) - fi* - (* 2 + \)/rUr = (m + *)*-, (5) 

where n is an integer called the meridional mode number. It determines 
the number of half periods of E in radial direction. The accuracy of 

(5) improves for large n, but is in most cases surprisingly good even for 
small values of n. Equation (5) permits an evaluation of the propaga- 
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Fig. 3 — Sketch defining regions of periodic and aperiodic field characteristics of a 
mode of azimuthal order v. 



tion constant fi for given mode numbers /x and v. This will now bo done 
for the parabolic ring structure sketched in Fig. 2. 



III. GROUP DELAY AND IMPULSE RESPONSE 

Figure 2 shows a cross-sectional view of the circular symmetric index 
distribution. The index has a maximum value n at r *= i? , decreases as 

n(r) = >?o[l - A(r - B ) 2 /« 2 ] for R - a < r < It + a, (6) 

and has a constant value 

n(r) = M (l - A) (7) 

everywhere else. We assume A to be small compared to unity, introduce 
the abbreviations 

fco = 2tt/i /X (8) 

and 

P =r - R , (9) 

and obtain, with the help of (1) and (6), 

fc»(r)«fcg(l - 2Ap 2 /a 2 ). (10) 
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In order to solve eq. (5) analytically, we assume in addition that 
Ro » a, which permits us to replace r by Ro in (5). As a result, 

(M + *)t = f Rl Vfcg - /3 2 - (. 2 + i)/i2g - 2A/cgpVa 2 dr, (11) 
which has the solution 

/3 = C^8 — (" 2 + \)/Rl - 2>/2A(m + |)Wa]». (12) 

The phase constant (3 of a propagating mode must furthermore 
fulfill the condition 

^ fco(l - A) (13) 

for the cladding field to have evanescent characteristics. This permits 
us to calculate the total number of propagating modes. Keeping v 
fixed, we first determine the number of modes m in a group with the 
same v. We do this by solving (12) for n with /3 = k (l — A). Since 
A« 1 

„ - ,,„„« + 1 = VA72 A - %±±£ + | ■ (14) 



This number decreases as v increases. The largest possible v is obtained 
for m = 1. Thus with (14) 

„ m « = (2fc8fl8A - *£! VA/2) 1 - i- (15) 

For the following approximations, we ignore the terms j. In this case, 
the sum over all m from v = to v max yields 

M = fafc^oA, (16) 

which is the total number of propagating modes. Using the same 
approximations, we can express m with the help of (15) in the form 

m = yfK/2ak (l - .7*0, (17) 

an expression which will be used later on. 

To calculate the mode delay, we first convince ourselves with the 
help of (14) and (15) that the n- and .-terms in (12) are small (of the 
order A) compared to kl. We therefore approximate |8 by 

p = k - (u 2 + i)/2fcoBg - V2A(m + \)/a. (18) 

The differentiation of /3 with respect to the radial frequency 
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Fig. 4 — Impulse response of the thin parabolic tube structure ; (a) individual mode 
groups, (6) power distribution for large mode numbers {Ro » a). 



co = cko/iio yields the group delay 



- 5* [i + (f« + £)/»»]. 



(19) 



where L is the fiber length and c the velocity of light in free space. 
Since t depends only on v but not on n, mode groups with the same v 
have the same delay. Consequently, if all modes are excited by equal 
pulses of unit energy at the fiber input, the output consists of pulses 
of energy m(v) delayed by t{v). If we ignore the delay Liio/c common to 
all modes and then insert (15) into (19), again neglecting the terms \, 
we can write the delay in the form 



, N , Imq Ln A v 2 

T(v) = t — = — 5- • 

C « ''max 



(20) 



Figure 4 illustrates the output distribution for the case in which the 
pulses are so narrow that individual groups are resolved. All pulses 
have the same (very small) width, and their heights correspond to the 
total energy m in each group. 

More meaningful than this plot is a plot of the energy per unit time 



p(r) = mdv/dr 



(21) 



which coincides with the power distribution in the case of very large 
mode numbers. We find dr/dv by differentiating (20) and, if the term 



IMPULSE RESPONSE OF FIBERS 



1167 



I.U 

in 






p 






z 






D _ 






„ 0.8 






> 






o: 






< 






a 






h 






E 0.6 


_ 




< 






z 






w 






g0.4 






a. 
05 




R ° \\\ 


111 




■ c"^..^^ 


cc 




a \J^ 


111 

w 0.2 






3 






a. 






| 











L J \ 1 ^sT" — --^ 



0.2 



0.4 



0.6 



0.8 



1.0 



RELATIVE DELAY rc/An n L 



Fig. 5 — Impulse response of the parabolic tube structure for various radius-to- 
thickness ratios. 



j in (14) is again ignored, we have 
caR ko 



_ caRpkl T / ALn y ( cr \*1 

This function is also plotted in Fig. 4. We find half the total output 
concentrated in the time interval 



T = 0.l2An L/c. 



(23) 



The remaining power is drawn out in a tail of length An L 'c. This tail 
is caused by modes of high azimuthal order v, modes which are not 
present in the 2-dimensional structure and are essentially equalized in 
the case of the concentric, parabolic profile (Fig. 1). In this respect, the 
parabolic ring structure is inferior to the corresponding concentric 
profile, yet an effective width of 0.12 An a L/c may be a useful improve- 
ment in comparison to a guiding structure with uniform index no which 
theoretically produces a width An u L/c. 

The condition R n » a was necessary for an analytic solution of the 
integral (11). Exact numerical results for arbitrary ratios R»a are 
shown in Fig. 5. The corrections with respect to (22) are largest for 
high azimuthal orders. An exact analytical solution can only be found 
for the maximum delay r{v max ) which becomes 
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As an example, consider a parabolic ring whose half width, a, is 
equal to the central radius Ro. Let n = 1.5 and assume A = 1 percent. 
The total width of the impulse response after 1 km of this fiber would 
be r(j'max) = 25 ns according to (22), but half the power is concentrated 
within the first 6 ns. Since high-order modes are usually lossier than 
the low orders, it is likely that much of the pulse tail does not reach 
the fiber end. 

IV. CONCLUSIONS 

The WKB approximation yields a simple characteristic equation for 
the propagating modes in fibers with arbitrary circular symmetric index 
distribution. We use this method to compute the impulse response of a 
fiber with a ring-shaped parabolic index distribution. We find that this 
structure has equalizing properties similar to the concentric parabolic 
index distribution, except for certain azimuthal mode orders, which 
lag behind, forming a rather long pulse tail. The rest of the power is 
concentrated in a time interval which, for a 1-km length and a relative 
index difference of 1 percent, is only 6 ns. 
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